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. A self-consistent analysis of K ^ 2tv and K ~* 3tv decays within a unique framework of 

' chiral dynamics applied to the QCD-corrected weak nonleptonic quark lagrangian has been 

performed. The results on K ^ 2tv amplitudes at O(p^), including the value for e /e, are 
compared with experiment to fix phenomenological B-factors for mesonic matrix elements 
' of nonpenguin and penguin four-quark operators. The dependence of these B-factors on 

^ , different theoretical uncertainties and experimental errors of various input parameters is 

investigated. Finally, we present our estimates at 0{p^) for the CP asymmetry of linear 
' slope parameters in the — > Stt Dalitz plot. 
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o 

' 1 Introduction 

^ il The starting point for most calculations of nonleptonic kaon decays is an efFec- 

1^ ■ tive weak lagrangian of the form E] 



Cl{\^S\ = l)^^/2GFVudV:,y^C.O., (1) 



> 

X 

5t 1 which can be derived with the help of the Wilson operator product expansion 

(OPE) from elementary quark processes, with additional hard gluon exchanges. In 
the framework of perturbative QCD the coefficients Ci are to be understood as 
scale and renormalization scheme dependent functions. There exist extensive next- 
to-leading order (NLO) calculations O U in the context of kaon decays, among 
others. These calculations are based on the possibility of factorization of short- and 
long-distance contributions into Wilson coefficient functions C^"^^ (/i) and mesonic 
matrix elements of four-quark operators Oi, respectively. The latter, however, can 
presently be obtained only by using nonperturbative, i.e. model-dependent, meth- 
ods, because not only perturbative QCD breaks down at scales fj, < IGeV, but 
also the QCD degrees of freedom (quarks and gluons) have to be replaced by the 
mesonic ones. 

Usually, the results of calculations are displayed with the help of i?-factors (bag 
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parameters) in the form 

vac. sat. ■) (2) 

i 

where the mesonic matrix elements of four-quark operators are approximated by 
their vacuum saturation values, which are real and /i- independent. In principle, 
factors should be estimated by some higher-order calculations in the long- 

distance regime, for instance, in l/A''c-expansion |S] in the form 1 -I- 0(l/iVc), or 
from the lattice approach. The preliminary stage of these calculations is best char- 
acterized by the long standing difficulties to explain quantitatively the well-known 
A/ = 1/2 rule. Recent lattice calculations seem to succeed in this respect, 
but are at the same time at variance (even in sign) with experimental values of 
Re(e je). Of course, the severe difficulties of all calculations of long-distance effects 
from "first principles" restricts the predictive power of (Q), leaving only the possi- 
bility of some semi-phenomenological treatment |31IZIIH1> which cannot be used to 
test the Standard Model (SM) in the kaon sector in a fundamental way. 

The main aim of the present paper is a further semi-phenomenological analysis 
of the long-distance (nonperturbative) aspects of the above lagrangian, especially 
in view of the actuality of the task to analyze the implications of the measured 
parameter of direct CP violation, e /e, on an alternative manifestation of direct 
CP violation by asymmetries in charged — > Stt decays. As the main result, 
we present numerical estimates of bag parameters and resulting asymmetries, in 
the form of probability densities, in order to demonstrate the experimental and 
theoretical uncertainties. 



2 General scheme of calculations 



In Oi are the four-quark operators, defined either by combinations of prod- 
ucts of quark currents (i = 1, 2, 3, 4, non-penguin diagrams) or, in case of gluonic 
{i = 5, 6) and electro-weak (i — 7, 8) penguin operators, by products of quark den- 
sities. The operators Oi with i — 1,2, 3, 5, 6 describe weak transitions with isospin 
change AI ~ 1/2 while the operator O4 corresponds to a A/ = 3/2 transition and 
operators Oy.g - to mixture of A J =1/2 and AI = 3/2 amplitudes. 

The operators Oi used here and in earlier work [Hl llOlITT] differ from those used 
in El El El and other papers H Oiniini (usually called Qi). Both sets are 
related by linear relations, which are given for easy reference below: 

Qi = 20i + JO2 + + , 
5 15 3 

Q2 - -2O1 + \02 + + , 

5 15 3 
Q3 = 2O1 + 2O2 , 

Q4 = -2O1 + 2O2 , 

Qs = 4a , = 2O5 + , 
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The general scheme of the meson matrix element calculation by using the weak 
lagrangian is based on the quark bosonization approach 10 . The bosonization 
procedure establishes a correspondence between quark and meson currents (densi- 
ties) and products of currents (densities). Finally, it leads to an effective lagrangian 
for nonleptonic kaon decays in terms of bosonized (meson) currents and densities: 

The meson currents/densities J^/fj^ and J£/^ are obtained from the quark 
determinant by variation over additional external sources associated with the cor- 
responding quark currents and densities jlOj . From the momentum expansion of the 
quark determinant to 0(p^") one can derive the strong lagrangian for mesons Ceff 
of the same order and the corresponding currents and densities J^/j^^ and J^^j^ to 
the order 0{p^"~^) and 0(p^"~^), respectively. Thus, the bosonization approach 
gives us the correspondence between power counting for the momentum expansion 
of the effective chiral lagrangian of strong meson interactions, 

and power counting for the meson currents and densities: 

^ jiP"'") (currents) ; C'f"'> ^ J^p"'^'> (densities) . 

Some interesting observations on the difference of the momentum behavior of 
penguin and non-penguin operators can be drawn from power-counting arguments. 
The leading contributions to the vector currents and scalar densities are of O(p^) 
and 0(p°), respectively. Since in our approach the non-penguin operators are con- 
structed out of the products of currents J£^, while the penguin operators are prod- 
ucts of densities J£, the lowest-order contributions of non-penguin and penguin 
operators are of O(p^) and 0(p°), respectively. However, due to the well-known 
cancellation of the contribution of the gluonic penguin operator at the lowest 
order ^H], the leading gluonic penguin as well as non-penguin contributions start 
from 0{p^) ^). Consequently, in order to derive the currents which contribute to 
the non-penguin transition operators at the leading order, it is sufficient to use the 
terms of the quark determinant to 0{p^) only. At the same time the terms of the 
quark determinant to 0{p'^) have to be kept for calculating the penguin contribu- 
tion at O(p^), since it arises from the combination of densities, which are of 0{p'^) 
and 0{p^), respectively. In this subtle way the difference in momentum behavior is 
revealed between matrix elements for these two types of weak transition operators; 
it manifests itself more drastically in higher-order lagrangians and currents. 

^ ) There is no cancellation of the contribution of the electromagnetic penguin operator Og at 
the lowest order and the leading contributions start in this case from 0{p°) 
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This fact makes penguins especially sensitive to higher order effects. In particu- 
lar, the difference in the momentum power counting behavior between penguin and 
non-penguin contributions to the isotopic amplitudes of if — > Stt decays, which 
appears in higher orders of chiral theory, leads to a dynamical enhancement of the 
charge asymmetry of the Dalitz-plot linear slope parameter [SI E| . 

In our approach the Wilson coefficients Ci in the effective weak lagrangian ^ 
are treated as phenomenological parameters which should be fixed from experiment. 
They are related with the Wilson coefficients Cf^'^ {^jl), calculated in perturbative 
QCD [7], via the S^-factors at the fixed renormalization scale /i = 1 GeV: 

The coefhcients Cf~^^ (/i) can be written in the form of a sum of z and y compo- 
nents, 

Cr^{,)^c[^\,)+rC\y\^.), r.-|^, (3) 

where only the y component contains the combination of CKM-matrix elements t, 
which, according to the SM, is responsible for CP violation - it has a non- vanishing 
imaginary part. For the amplitudes of nonleptonic kaon decays follows the same 
decomposition into z- and y- components, 

Ai - + tA^^ , (4) 

where Ai are isotopic amplitudes oi K ^ 2-k decay corresponding to the tttt final 
states with the isospin / — 0,2. 

With the dominating contributions A^p from the four-quark operators O^, the 
K — > (27r)/ amplitudes may be written as 

In the calculation of the A^j \ sizable corrections connected to isospin symmetry 
breaking have been taken into account (see below). 

The observable effects of direct CP violation in the nonleptonic kaon decays are 
caused by the y-components of their amplitudes (@J). In particular, the ratio e / e 
can be expressed as 

Ref-^ =ImAt(Po-i'2), Pi ^ ^ 4-v > 

where ImA* = \^V:,Vtd = |Kb| |T46|sin,5; w = 4'V4"^- 
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3 Theoretical and phenomenological uncertainties 

There are the following theoretical and phenomenological uncertainties which 
appear both from short distance (Wilson coefficients) and long distance (effective 
chiral lagrangians and _Bi-factors) contributions to the kaon decay amplitude: 

— Dependence on the parameter Im r ^ Im At arising fr'om the imaginary part 
of Wilson coefficient © . 

— Regularization scheme dependence which arises when calculating Wilson co- 
efficients beyond the leading order (LO) of QCD in the next-to-leading or- 
ders: naive dimensional regularization (NDR), t'Hooft-Veltman regularization 
(HV). 

— Dependence of Wilson coefficients on choice of the renornialization point ^, 
taken below as 1 GeV, and QCD scale Ap^, contained in the interval (325 ± 
110) MeV. 

— Factors Bi [i — 1, 4, 5, 8) for dominating contributions of four-quark operators 
Oi to the meson matrix element (jSJ- 

— Dependence of meson matrix elements on the structure constants L2, is, £4, 
-L5, of the general form of the effective chiral lagrangian introduced at 
0{p^) by Gasser and Leutwyler |l(ij . 

— Dependence of meson matrix elements on the structure constants of the ef- 
fective chiral lagrangian at O(p^) [I7j . 

— Dependence on choice of the regularization scheme to fix the UV divergences 
resulting from meson loops. 

In the present paper we have combined a new systematic calculation of mesonic 
matrix elements for nonleptonic kaon decays from the effective chiral lagrangian 
approach with Wilson coefficients Cf^^{^), derived by the Munich group |Zj for 
/Lt = 1 GeV and mt — 170 GeV. For the parameter ImAt we have used the result 
obtained in ^H]: ImAt = (1.33 ± 0.14) • lO"'*. We performed a complete calcu- 
lation of if ^ 27r and K in amplitudes at 0{p^) including the tree level, 
one- and two-loop diagrams. For the structure coefficients Li of the effective chiral 
lagrangian at 0{p'^) we used the values fixed in jl9| from the phenomenological 
analysis of low-energy meson processes. The structure coefficients of the effective 
chiral lagrangian at O(p^) have been fixed theoretically from the modulus of the 
logarithm of the quark determinant of the NJL-type model (see |21 for more de- 
tails). The superpropagator regularization has been applied to fix UV divergences 
in meson loops, thereby the renormalization scale was /i5p = 1 GeV. The isotopic- 
symmetry-breaking corrections which arises from both the ir^-r] and i^^-r] mixing 
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and the final-state tt^ — 7r° mass difference were taken into account ^) . In one 
can find more technical details of the calculation of K ^ 2tt amplitudes. 

4 Phenomenological analysis o{ K 2% decays 

In our phenomenological analysis the results on the K ^ 2tt amplitudes, includ- 
ing the value for s /e, are compared with experiment to fix the phenomenological 
i3i-factors for the mesonic matrix elements of nonpenguin and penguin four-quark 
operators. As experimental input we used the experimental values of the isotopic 
amplitudes ^0*^2^'* fixed from the widths of K ^ 2tt decays and the world average 
value Res' /e = (16.2 ±1.7) x 10"^ which includes both old results of NA31 HI] 
and E731 [211 experiments and recent results from KTeV 23^ and NA48 The 
output parameters of the performed K ^ 2tt analysis are the factors Bi, B4 and 
B5 for a fixed value of Bs- 

The dependence of the i?i-factors on different theoretical uncertainties and 
experimental errors of various input parameters is investigated by applying the 
"Gaussian" method. Using Wilson coefhcients derived in |Zj in various regular- 
ization schemes (LO, NDR, HV) for different values of the QCD scale A.^1^ ,we 

calculated the probability density distributions for Sj-factors obtained by using 
Gaussian distributions for all input parameters with their errors. As an example, 
the probability densities for the parameters Bi, B^^ B^ calculated with = 1 
and = 325 MeV are shown in Figs. ^ [3 and 01 Upper and lower bounds for 

Bi-factors (i = 1,4,5) for different values of a!^^ in LO, NDR and HV regular- 

ization schemes (Bg, = 1) obtained by the Gaussian method are shown in tabled 
The limits without parentheses correspond to the confidence level of 68% while the 
limits in parentheses - to the confidence level of 95%. 

Fig. 131 shows the necessity for a rather large gluonic penguin contribution to 
describe the recently confirmed significant experimental e value (the factor B^ is 
found well above 1). It should be emphasized that the dependence of Bi (i — 1, 4, 5) 
on i?8 is very small even within a wide range of its values < Bg < 10. Therefore 
even for Bg = values of -65 > 2 are necessary to explain the large value of e /e. 

Even for larger values of B^, the contributions of nonpenguin operators to the 
AI —1/2 amplitude are still dominating (see Fig. ^ and also upper and lower 
bounds for the relative contribution of penguin operators to the AI =1/2 ampli- 
tude in table EJl. The large Bi and B^ values may be a hint that the long-distance 
contributions, especially to AI = 1/2 amplitudes, are still not completely under- 
stood. An analogous conclusion has been drawn in ^HI> where also possible effects 
from physics beyond the Standard Model are discussed. 

^) In 111! the ■K'^-rj-'q mixing contribution to the isospin breaking has been computed only at 
tree level. The importance of the correction due to the final-state tt^ — vr'' mass difference was 
discussed for the first time recently in I2UI . This mass difference leads to a sizable correction to 
e /e, in spite of its smallness. 
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B5 



Fig. 3. Probability density distributions for factors with Bg = 1 



Czech. J. Phys. 53 (2003) 



A9 



A. A. Belkov, G.Bohm, and A.V. Lanyov 



Table 1. Upper and lower bounds for Bi factors (i = 1, 4, 5). The limits without paren- 
theses correspond to the confidence level of 68% while the limits in parentheses - to the 

confidence level of 95%. 



B, 


AllL 

MS 

(McV) 


LO 


NDR 


HV 


min max 


min max 


min max 


Bi 


215 


5.4 6.4 

( 5.0 6.9 ) 


6.1 7.2 

( 5.7 7.9 ) 


5.8 6.8 

( 5.3 7.4 ) 


325 


4.9 5.8 
( 4.6 6.3 ) 


5.7 6.8 
( 5.3 7.4 ) 


5.2 6.2 
( 4.7 6.7 ) 


435 


4.5 5.3 
( 4.2 5.8 ) 


5.3 6.3 
( 4.9 6.8 ) 


4.6 5.4 
( 4.3 5.9 ) 


Bi 


215 


0.68 0.71 

( 0.67 0.72 ) 


0.61 0.63 

( 0.60 0.65 ) 


0.64 0.66 

( 0.62 0.68 ) 


325 


0.72 0.75 
( 0.70 0.76 ) 


0.64 0.66 
( 0.63 0.67 ) 


0.67 0.70 
( 0.66 0.71 ) 


435 


0.76 0.79 
( 0.75 0.81 ) 


0.66 0.68 
( 0.65 0.70 ) 


0.73 0.76 
( 0.72 0.77 ) 




215 


2.5 3.9 
( 2.0 4.9 ) 


2.8 4.5 
( 2.3 5.6 ) 


3.5 5.3 

( 2.9 6.7 ) 


325 


1.8 2.9 
( 1.5 3.7 ) 


2.0 3.2 
( 1.6 4.0 ) 


2.3 3.6 
( 1.8 4.5 ) 


435 


1.4 2.2 
( 1.1 2.9 ) 


1.3 2.0 
( 1.0 2.7 ) 


1.5 2.5 
( 1.2 3.1 ) 



5 Direct CP violation in K ^ 'in decays 

Finally, the predictions for the CP asymmetry of linear slope parameters in the 
— > Stt Dalitz plot are discussed. These predictions are based on a new calcula- 
tion oi K ^ "in amplitudes at 0(p^) within the same cfFcctivc lagrangian approach. 
The obtained if ^ Stt amplitudes include the jame theoretical uncertainties as in 
case of the K ^ 2-k analysis. The values of i?4, fixed from the — > 27r 
analysis are used as phenomenological input to the if ^ Stt estimates which have 
been performed in a self-consistent way by the Gaussian method. 

The linear slope parameter g of the Dalitz plot for if — > Stt decays is defined 
through the expansion of the decay probability in the kinematic variables 

|T(if ^ 37r)|^ cx l+gY+ ... 

where F is a Dalitz variable: 

^ sz- So 2 rn^K , 2 
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aJ^s = 325 MeV 






LO 






NDR 






HV 


1 1 1 .. -r' 1 


IV 1 1 . ^ u - 'I' 1 





0.02 0.04 0.06 0.08 0.1 0.12 

A„(p)/Ao(tot) 

Fig. 4. Probability density distributions for the relative contribution of penguin operators 

to the A J ~ 1/2 amplitude 



Table 2. Upper and lower bounds for the relative contribution of penguin operators to 
the AI = 1/2 amplitude. The limits without parentheses correspond to the confidence 
level of 68% while the limits in parentheses - to the confidence level of 95%. 





MS 

(MeV) 


LO 


NDR 


HV 


min max 


min max 


min max 




215 


0.023 0.029 
(0.019 0.031) 


0.062 0.078 
(0.053 0.084) 


0.031 0.039 
(0.027 0.042) 


325 


0.025 0.031 
(0.021 0.034) 


0.071 0.088 
(0.061 0.096) 


0.032 0.040 
(0.027 0.044) 


435 


0.026 0.033 
(0.022 0.036) 


0.079 0.098 
(0.067 0.107) 


0.043 0.054 
(0.037 0.059) 



and the index tt^ belongs to the odd pion in the decays n^Tr^n^ and 

TT^TT^Tr*. Direct CP violation leads to a charge asymmetry of the linear 
slope parameter, 

Fig.Elshows the probability density distributions for tt^tt^tt^ and 

tt^tt^tt'^ decays calculated with _Bg = 1 and A^IL = 325 MeV. Upper and low 

bounds for Ag^^_ and Agoo+ for different values of A^IL in LO, NDR and HV 

AI S 

regularization schemes (_Bg = 1) obtained by the Gaussian method are shown in 
tableOl Summarizing these results, we have obtained the following upper and lower 
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Fig. 5. Probability density distributions for the CP asymmetry of linear slope parameters 
of — » TT^TT^TT^ and — > 7r°7r°7r* decays 



bounds for the charge symmetries of the Hnear slope parameter: 

1.6 < Ag++_ ■ 10^ < 4.2, 0.9 < A500+ • lO'^ < 2.2 with CL=68%; 

1.1 < A5++_ • W < 7.6, 0.6 < A500+ • W < 3.3 with CL=95%. 



When comparing these results with those found above for the phenomenological 
-B-factors, one finds a much reduced dependence on the regularization scheme and 
on the scale aIIL. Once more this type of asymmetry ratios turns out to be more 
stable not only against systematic experimental errors (efficiencies, among others), 
but also with respect to theoretical uncertainties of its parts. 
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Table 3. Upper and low bounds for Ag++- and /S.goo+ (in units 10~*). The limits without 
parentheses correspond to the confidence level of 68% while the limits in parentheses - to 

the confidence level of 95%. 





MS 


LO 


NDR 


HV 




(MeV) 


min 


max 


min 


max 


min 


max 




2f5 


f.6 


4.f 


f.5 


3.6 


f.6 


4.f 






( 1-1 


7.6 ) 


(f.O 


6.0 ) 


( 1-1 


7.4 ) 


A(?++- 


325 


f.6 


4.2 


f.5 


3.6 


f.6 


4.1 






( 1-1 


7.6 ) 


(f.O 


5.9 ) 


( 1-1 


7.3 ) 




435 


f.6 


4.2 


f.5 


3.5 


f.6 


4.0 






( 1-1 


7.6 ) 


(f.O 


5.7) 


(f.f 


7.0 ) 




2f5 


0.9 


2.2 


0.9 


2.f 


0.9 


2.2 






( 0.6 


3.3 ) 


( 0.6 


3.f ) 


( 0.6 


3.3 ) 


Agoo+ 


325 


0.9 


2.2 


0.9 


2.f 


0.9 


2.2 






( 0.6 


3.4 ) 


( 0.6 


3.f ) 


( 0.6 


3.3 ) 




435 


0.9 


2.2 


0.9 


2.f 


0.9 


2.2 






( 0.6 


3.4 ) 


( 0.6 


3.0 ) 


( 0.6 


3.3 ) 



6 Conclusion 

With some experimental updates and theoretical refinements our new estimates 
confirm the dynamical enhancement mechanism for the charge asymmetry by 
higher order contributions in the effective chiral lagrangian approach which was first 
observed in [S] . The predicted slope parameter asymmetry, although small, may be 
in reach of current high statistics experiments [53]. Already with lower statistics, 
new measurements of quadratic slope parameters oi K —f 3tt decays, including 
the neutral channels, would lead to an improved theoretical understanding of the 
nonperturbative part of nonleptonic kaon decay dynamics. 

For £ /e, the results of analysis are very sensitive to final-state interactions, 
isotopic-symmetry-breaking effects and other refinements of the model calculations. 
Despite of considerable theoretical efforts dedicated to the calculation of Wilson co- 
efficients and bag factors, the remaining uncertainties are still very large (compare 
also |2fij . where a similar conclusion is reached, when taking into account the full 
range of all input parameters and all theoretical uncertainties). The fact that phe- 
nomenological values of -B^-factors {i — 1, 4, 5) considerably differ from unity shows 
that the long-distance contributions are still not completely understood. Besides 
the (rather remote) possibility, that the SM has to be revised in the kaon sector, 
there are further, more exotic developments in long-distance QCD going on |27|. 
which may have manifestations also here. Another very active direction is the inves- 
tigation of the problem of matching between short- and long-range renormalzation 
schemes by means of (effective) color singlet boson models. For a recent study. 
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especially of CP violation in K Stt decays, see [5S| and papers cited there. 
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